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Abstract. We calculate perturbatively the multifractality spectrum of wave- 
functions in critical random matrix ensembles in the regime of weak 
multifractality. We show that in the leading order the spectrum is universal, 
while the higher order corrections are model-specific. Explicit results for the 
anomalous dimensions are derived in the power-law and ultrametric random 
matrix ensembles. 



PACS numbers: 72.15.Rn, 71.30.+h, 05.45.Df 



Universal and non-universal features of the multifractality exponents. 



2 



1. Introduction 



Critical wave- functions tp(r) at the point of the Anderson localization transition [T] 
are usualljQ multifractal [3] . Their statistical properties can be characterized by a set 
of exponents d q , which describe the scaling of the wave- function moments with the 
system size L in the limit L — > oo. Using the bar to denote the disorder averaging, 



one can define d q via the relations 



I q = d d r |V(r)| 2 <2 oc L"^^" 1 ) . (1) 



For completely extended states d q = d; for completely localized states d q — 0. The 
hallmark of multifractality is a g-dependent d q , different from both and d. 

Although detailed characterization of the critical wave-function multifractality 
remains an active area of research (see [T] ) , relatively few explicit analytical results for 
d q are available in the literature. One such result is due to Wegner [4], who calculated 
the anomalous multifractal dimensions A g = {d q — d)(q— 1) for the Anderson transition 
in 2 + e dimensions: 

A 9 = q(l - q)e + ^q(q -l)(q 2 -q + l)e 4 + 0(e 5 ), /3 = 1, 

A q = q (l-q)^) 1/2 - 3 -^q\q-lfe* + 0(e^), = 2, (2) 

with standing for the Riemann zeta-function. These expansions were derived in 
the replicated nonlinear er-model with a coupling constant t oc e in the orthogonal 
symmetry class (time-reversal invariant systems, (3 = 1) and t oc y/e in the unitary 
symmetry class (systems with broken time- reversal invariance, /3 = 2). Note that the 
second and third order corrections are equal to zero in both cases. 

Another explicit result was derived more recently [3 [6] in the framework of 
power-law banded random matrices (PLBRM), introduced in [7] as a useful toy-model 
with multifractal eigenstates. In the regime of weak multifractality the anomalous 
exponents are given in the lowest order bj|§] 

A q = q(l-q)t + (f), (3) 

where t//3 plays the role of the coupling constant in the effective a- model and is 
assumed to be small. Comparing eq. ([2]) with eq. ([3]) we notice that the leading 
orders are the same if one identifies e ( \J ej2 for [3 = 2) with the coupling constant t//3 
of the corresponding cr-model. This observation can be generalized to an arbitrary 
model at criticality in the regime of weak multifractality (see Section [2]) , making the 
leading order contribution in A q model-independent or universal. While the critical 
exponents of short-ranged models are numbers, in the long-ranged models they are 
functions of parameters (e.g. of bandwidth in PLBRM). Universality then means that 
the functional forms of critical exponents from different models can be mapped onto 
each other. Remarkably, in PLBRM, a similar universality also exists in the opposite 
regime of strong multifractality [8]. 

The aforementioned universality poses a question: does the knowledge of higher 
order contributions to A q allow to distinguish among different models? These 
contributions may contain model-specific information, making A q different for different 

| For a rare exception to this rule see [2]. 

§ In this work t is 87r times smaller than in the cited reference. 
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models within the same symmetry class. One should stress, however, that since A q is 
a critical exponent, it is expected to be insensitive to any short-scale deformation of 
a model within its universality class. 

In the present work we calculate the anomalous multifractal dimensions to the 
second order in the coupling constant for two critical random matrix ensembles: 
PLBRM and the ultrametric random matrices [8]. In the unitary symmetry class 
(0 = 2), we find that the second order correction is equal to zero. It follows from our 
derivation that the absence of the second correction is very general and not restricted 
to the particular models we consider. In the orthogonal case (/? = 1), the second 
order correction turns out to be non-zero and different for the two models considered. 
This fact should be contrasted with the situation in 2 + e dimensions and may be a 
consequence of the long-range nature of the model Hamiltonian. In addition, we find 
that the first deviation of the multifractal spectrum from its parabolic shape occurs 
in the third order correction and is also model-specific. 

2. Perturbation theory in the non-local non-linear cr-model 

Non-linear supersymmetric cr-model is the standard field-theoretical description of 
a single quantum particle in a disordered medium [9]. When computing critical 
exponents, insensitive to short-scale details, one can take the equivalent cr-model 
description as the starting point instead of using the underlying random matrix- 
type Hamiltonian explicitly. From this angle, let us consider a ID chain of L 
sites enumerated by a discrete spatial coordinate r, and associate with every site 
a supermatrix Q r . The effective cr-model is defined by its action: 

S[Q] = ~^- t Yl MQrQv) -ivYl Str (a-A). (4) 

r,r' r 

Here the kernel U r y defines the non-local term in the action ( "kinetic energy" ) , and 
can be viewed as a symmetric L x L matrix with a zero-mode: ^ r U r ^ r i — 0. The 
coupling constant £ 1 is the small parameter of the perturbative expansion, and (3 is 
the usual symmetry parameter of the random matrix theory: /3 = 1 in the orthogonal 
symmetry class (Q being a 8 x 8 matrix) and /3 = 2 in the unitary symmetry class (Q 
being 4 x 4). The matrix A = diag(l,— 1) in the advanced-retarded representation, 
and n > represents the effective level broadening. 

Let us sketch the main steps of the standard perturbative treatment of the model 
(see PHIE] and references therein for more details). The wave- function moments (fTJ) 
(also known as the inverse participation ratios or IPR's) can be extracted as 

I q = lim^/ VQe™, e F ^ = (V st r«(Q r Ak)) e~ s , (5) 

J r 

where k = 1 ® diag(l,— 1) in the advanced-retarded representation. Following |10j . 
we will use the rational parametrization of the matrix Q(r) with the zero-mode Qq 
separated: 

Qr = TAg(W r )T-\ Q Q = TAT -1 , g(W) = (1 + ±W)(1 - \w)^ . (6) 

Here T is a position-independent matrix and the spatial variations are represented by 
W r - a small ( 0{yi)) block-off-diagonal matrix with a zero mode: ^ r W r = 0. The 
change of variables Q — > {Qo, W} in the functional integral incurs a Jacobian, known 
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to be J = exp ^ J2 r W 2 + ^(i 2 ) in the orthogonal case, and 0(t 2 ) in the unitary 

case. The perturbation series is based on the Taylor expansion g(x) = 1 + x + ^.x 2 + ... 
Substituting such an expansion into the action one separates the O(l) part of the 
"kinetic energy": So = ^ J2rr> U r y str(W r W r ') and treats all the other terms in ([5]) 
as a perturbation. The averages (• • •) with respect to S a are computed by the Wick 
theorem based on the elementary contractions: 

£ = 2: ( str(W r P) str(W r /P)) = tU ry str{PR - PARA), 

(str(W r PW r /P)) = tfI r>r <(strPstrP- str(PA) str(PA)), (7) 



13 = 1: ( str(W r P) str(W r /P)) = flT r , r * str(PP - PARA + PARA - PR), 
(str(W r PW r 'P)) = ffl r , r <(strPstrP- str(PA) str(PA) + str(PAPA - PR) ) . (8) 



gfefdQo] 



Here P and i? are arbitrary matrices and R — M T R T M, where M is an orthogonal 
matrix of the time-reversal transformation. The propagator n r r /, associated with the 
action So, is defined by the identity C/ rijr n r . jr2 = <5 rij r 2 — ^, where the subtraction 
is required by the zero-mode of U r ^ r > ; it also ensures that the propagator has a zero- 
mode as well: n rl .< = 0. Moreover, we consider spatially uniform systems in 
which the diagonal elements U r _ r and n r r are independent of r. Implementing such a 
perturbation theory, one effectively integrates out all the factors W r containing non- 
zero momenta, and is left with the integral over the zero-mode Qo only, which is 
simplified in the limit rj — > and can be done using any standard parametrization 
of this matrix 9 . As a result, one can prove in the unitary case the following 
correspondence: given that 

Lstr 9 (Q Afc) A, + Ai(eLstr(Q A)) +... + A n (eL str(Q A)) " , (9) 

where Ai are some constants (in the present calculation only the first three are non- 
zero), the zero mode integrals over Qq yield 

I q = ^(Ao + Ma-l) + ...+A ^ + n _ 2) f- ), (10) 

and a similar expression with the prefactor (2q — 1)!!/L 9_1 in the orthogonal case. 

Performing a systematic expansion up to the order t 2 in this way, we find for the 
unitary case: 

IP - #r (i+^n^+^g^nL-t 2 ^- 1 ^- 1 ^ ,)^^), (ii) 

where H 2 = L~ 2 J2rr' n 2 r /. To obtain this result, one needs to keep track of terms 
with up to four WV-factors everywhere except in the expansion of the "kinetic energy" 
in S, where one has to go up to six W r 's because of the 1/t prefactor. Remarkably, in 
the unitary case these extra terms vanish after averaging over W r , as does the 0(t 2 ) 
contribution of the Jacobian J. The latter statement is still true in the orthogonal case, 
but the former is not: due to a more complicated structure of the W-W contractions 
in this "kinetic" contribution K survives the average: 

r = (i + Ms-DiW + * 2 ^^n* r + t 2 fciln 2 r _ 

t2 g ( g -i)(4 g - 5 ) n2 + _ \ + 
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K = L- 1 ^ U r yTl r yIl r y,Il r ,y, -Il 2 rr ~(L-l)Il 2 . (13) 

r,r' ,r" 

Also, unlike the unitary case, the Jacobian J has an 0{t) part that contributes to the 
penultimate term in the brackets of (|12l) . 

In the limit of weak multifractality, the anomalous dimensions A q are 

parametrically small and hence I q oc L 1 ~ q exp(— A 9 logi) = L 1 ^ q (^l — A q \ogL + 

^A 2 log 2 L + 0(A q )). Having found the large L asymptotic behavior of equations (|TT|) 
and (I12p , we can check that it is indeed of the required form and extract A q up to the 
order t 2 . Note that in order to demonstrate the first proposition it is sufficient to show 
that for large L (i) the kinetic contribution behaves as K — > —ill 2 ,, + ci logL + c 2 , 
and (ii) II2 tends to a constant. The latter appears to be a universal fact: if II2 grew 
with L, it would contribute to the anomalous dimension. However, the prefactor in 
front of H2 does not have the symmetry q — > 1 — q, which the anomalous dimension 
of any a- model should possess on general grounds [12]. Therefore, we can see already 
from eq. (|11|) that in the unitary case the second order correction to the anomalous 
dimension vanishes. Below we will apply the results of this section to two types of 
critical systems: the a-model corresponding to the critical PLBRM ensemble [7] as 
well as the critical er-model, which can be derived from the n-orbital version of the 
ensemble of ultrametric random matrices introduced in [8] . 



3. Critical Power-law Banded Random Matrix ensemble 

Let us briefly recall the definition of the critical PLBRM and the corresponding a- 
modcl. Consider a system of L sites whose Hamiltonian is an L x L random matrix 
r—r'\), and further confine it to a circle by imposing periodic boundary 
conditions. The random matrix G is taken from the standard Gaussian ensemble 
(unitary or orthogonal), and a(r) is a function that controls the spatial decay of the 
disorder: a(r) ~ 1 for r < b and a(r) — b/r asr-> 00. In the limit of large bandwidth 
b ^> 1 this model can be mapped [7] onto the non-local cr-model Q with 

U r , r > = (A~ l ) r ,r< ~ 6 r ,r> Aj) » * = < 1, (14) 

i 

where we introduced an L x L matrix A with entries A r>r i = a 2 (\r — r'\). The purpose 
of subtracting the diagonal term in U r _ r ' is to ensure ^ r V r y = 0. Translational 
invariance suggests working in the momentum space. For small momenta k <^ L, the 
discrete nature of the underlying lattice becomes irrelevant, and we can safely employ 
the long- wave asymptotic of the Fourier components: Uk ~ |fc| and life ~ l/|fe|. In 
fact, since the anomalous dimensions must be insensitive to short-scale details of the 
model, only such small momenta are important and one can use the above asymptotic 
expressions as exact definitions for an effective model. Defining the cutoff fc m via 
L = 2k m + 1, we then arrive at: 

U rr ,=2L- 2 V e^^-^lkl, Tl rr , =\ V e J T^-'-')J_. (15) 
fe=±l,... 2 fe=±i,... 

In particular, H r ^ = J2k=i = l°g(^me 7 ) + 0(l/fc m ), where 7 « 0.577 is the Euler- 
Mascheroni constant, and H2 = 2"Sfc=i V^ 2 = fo + 0(l/k m ). These formulas are 



Universal and non-universal features of the multifractality exponents. 



G 



sufficient for our purpose in the unitary case: 

IP = #T (l + log L + eibzlL log* L + fc) + 0(t% (16) 

where only the leading terms in L — > oo are kept and by c we hereafter denote a 
constant independent of both L and t. The above expression implies 

\ = ^Y^t + 0(t 3 ), /3 = 2. (17) 

The second order correction is zero, as announced. The j3 = 1 case is less trivial, since 
eq. (|T3| also requires the asymptotic 

u r , r ,n r , r ,n r , r „n r ,, r „ = ^(y fc m + \n 2 r r - n r , r + c ) + ofc 1 ), (is) 

r.r' ,r" 

which is extracted using the Euler-Maclaurin formula [13] and the asymptotic 

2 2 

expansion of the digamma function at infinity. Note that the leading term ^-fc m ~ 
(which would be incompatible with the assumption of multifractality) cancels out with 
the last term in eq. (JT3J - As a result, 

I( q 0) = (l + *(1 - *)8(8 - 1) l°g^ + t 2 g2(g 2 ' 1)2 log 2 L + t\) + 0(f), (19) 

implying the anomalous dimension 

A q = q(l-q)(t-t 2 ) + 0(t% /3 = 1. (20) 

In view of the long-scale nature of the calculation, one can try to repeat it replacing 
all momentum sums with integrals. The values of individual computed quantities 
will change. For instance, one will have II2 = 1, rather than 7r 2 /12, which could 
conceivably destroy the vital cancelation of the 0(k m ) term in (|T3|) . However, this 
does not happen: the prefactor in front of k m in the first term of (fTSj) also changes 
so as to preserve the cancelation, and in this way the anomalous dimensions (|17|) and 
(|20[) are found correctly. This is a clear manifestation of the fact that multifractality 
is insensitive to the discreteness of the model. 



4. Critical ultrametric model on a binary tree 

The critical ultrametric random matrix model (with binary splitting at each level of the 
hierarchy) was introduced in [S] . The model itself can not be exactly mapped onto the 
nonlinear a-model. However, by associating with every site of the underlying lattice 
n internal degrees of freedom ("orbitals") and coupling them by Gaussian random 
matrices in a way respecting the original lattice structure (see e.g. |14j for a more 
detailed description of the procedure), one can derive an effective nonlinear a— model 
in the limit n — > 00. This a- model is expected to share the critical properties with the 
original microscopic model. Performing the necessary steps in the ultrametric case, we 
once again arrive at the non-local cr-model (j4|) with the kernel U r y expressed in terms 
of the variance matrix of the ultrametric random Hamiltonian A T y — (|-ff r .r'| 2 ) as 

U r . r ' = (A~ 1 ) r . r / — 8 r y (^2i Ai,jj . Now the system size can only be L = 2 k , where 
the integer k is the number of levels in the hierarchical construction. The variance 
matrix A is constructed recursively: 

^ ^(rtt!)' (21 » 
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where E is a square matrix with all entries equal to 1. 

Unlike the power-law case, the value of the coupling constant f in Q is not 
pre-determined by microscopic parameters, but can be chosen ad hoc. Being mainly 
interested in developing the perturbative treatment of the model, we consider this 
coupling to be small. Anticipating the result, it is actually convenient to replace t in 
(j4| with the combination 

~ = _61og2_ { ) 

In order to diagonalize matrices with ultrametric block structure, it is technically 
convenient to use the so-called replica Fourier transform (RFT), as presented in |15j . 
The quantities of interest are then written as discrete sums, which can be subsequently 
expressed in terms of the g-digamma function with q — 2, whose asymptotics are 
known [13]. For instance, we obtain II r , r = -(w + 2 - 2 1 ~ k ) £* =1 2 j ~ 1 " + lzl%~ k , 
whence Ti r ,r = ^7^2 + Proceeding in this way, we recover in the unitary 

case the same formulas (fT6|) and (|T7|) as in the PLBRM model. In the orthogonal case, 
we again need to take care of the contribution from the "kinetic" terms (fTB")) . After a 
rather tedious calculation we arrive at the final result: 

(?|zI2( 1+ ( t _ t ^.-,w 

A« = 9(l-?)(i-f 1: f log2)+0(t 3 ), 11 = 1. (24) 

In contrast with the power-law case, the anomalous dimensions for the ultrametric 
case could not be obtained by simply replacing sums with integrals. 

5. Non-parabolic corrections to multifractal spectrum in the unitary case 

We have found that in the first two orders the anomalous dimensions A g are given by a 
quadratic polynomial in q. This means that the multifractal spectrum /(a), obtained 
from d q by a Legendre transform [I], is parabolic. The first correction to A q which 
induces non-parabolicity of f(a) occurs in the next order 0(t 3 ), where we restrict 
our calculation to the unitary symmetry class (/3 = 2) for simplicity. Continuing the 
perturbation series, we find in the PLBRM case 



^ = ^7^fl + (^^^l°g2)^-l)logi+t 2 ^V L l°g 2 i+t 2 C )+0(i 3 ),(23) 



1 Li 



1 ' 3(? -i)„ , ^te-i) 2 ^ , M 3 (q-i)* 3 



1 + t — - ' Ur, r + ^ g Kr + * ~ 4 g Kr + 



t d A(L) q 2 (q - l) 2 + t A B(L) q(q - 1) + t A cj + 0{t 4 ). (25) 

In the ultrametric case, the expression is the same up to a change t t. A 
prefactor C = 1 — t 2 q ( q ~ 1 '>( 2q ~ 1 ') n 2 has been isolated. One can see that non- 
parabolic corrections to A q arise only from the term denoted A(L), which reads 

A{L) = -|n?, r + iz En.ni n ^ 2 + sT ^r u r 2 ,r 3 Ur 1 ,rjQ,, ri n* a , n - We computed the 
sums numerically in both models and fitted the result by a polynomial of logL. The 
absence of log 3 L and log 2 L terms can be stated with confidence: in our numerical 
calculations the coefficients in front of these terms were < 1CP 4 . These terms are 
expected to be absent, since their presence would be incompatible with the multifractal 
scaling. As for the linear term, it was found to be different in the two models: 

ft \ 3 

Yplbrm (£) = -J/plbrm log L + C, A ultra (i) = - (jj 2/ultra log L + C, (26) 
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where t/t was given in eq. (I22[) and 

J/ P ibrm = 0.407 ±0.004, j/uitra = 0.486 ±0.008. (27) 

The dependence of A(L) on a; in eq. (|26[) in the ultrametric case is natural: it means 
that A q depends on u> only via t. We also checked this fact numerically for three values: 
u> = 0, 0.5 and 1. As a result, the non-parabolic part of the anomalous dimension is 
given, both in PLBRM and ultrametric cases, by 

A™ n ~ pal = t 3 yq 2 (l-q) 2 + 0{t 4 ), (3 = 2, (28) 

with the value of y taken from (|2"T)) . The fact that this value is different for the 
two models, while in the first two orders their dimensions were identical (eq. (|17[) ). 
demonstrates the announced "non-universality" of the anomalous exponents beyond 
the parabolic approximation. 

Conclusion 

The second order corrections to the anomalous dimension for f3 — 1 (eqq. (1241) and 
(|20I0 appear different in the two models: it is roughly 50% larger in the ultrametric 
case (i&log2 ss 1.5018). A somewhat similar situation (identical first order, but 
different second order corrections) arises in the comparison of the classical Ising model 
with a power-law interaction and the hierarchical Dyson mode|jj] [16) . However, since 
the q-dependence of our A q is parabolic in the first two orders, this difference could 
be in principle absorbed by a suitable redefinition of the coupling constant in the 
ultrametric case. To resolve the question of alleged non-equivalence of the two models, 
and hence of the non-universality of the anomalous dimensions in long-ranged systems, 
we extended our analysis of the unitary case to the next order of the perturbation series 
and calculated the first non-parabolic contribution to A q . We observed a difference 
between the two models (eq. (|27p ). which can no longer be removed by a simple 
redefinition of couplings. 

Instead of a simple scaling law ([T]) , I q can in principle behave as a sum of several 
different powers of L, the highest of which dominates in the limit L — > oo and gives the 
multifractal spectrum. Nevertheless, up to the second order in the coupling constant 
our results are consistent with the single-power scaling (JTJ). This is true for any non- 
local cr-model and not only for the two examples considered here. 

IR and AO acknowledge support from the Engineering and Physical Sciences 
Research Council [grant number EP/G055769/1]. 
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